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Using data from the recent BONuS experiment at Jefferson Lab, which utilized a novel spectator 
tagging technique to extract the inclusive electron-free neutron scattering cross section, we obtain 
the first direct observation of quark-hadron duality in the neutron F 2 structure function. The data 
are used to reconstruct the lowest few (N = 2, 4 and 6) moments of F 2 in the three prominent 
nucleon resonance regions, as well as the moments integrated over the entire resonance region. 
Comparison with moments computed from global parametrizations of parton distribution functions 
suggest that quark-hadron duality holds locally for the neutron in the second and third resonance 
regions down to Q 2 « 1 GeV 2 , with violations possibly up to 20% observed in the first resonance 
region. 


I. INTRODUCTION 

Inclusive lepton scattering has for many decades been 
the most important tool for probing the internal quark 
and gluon (or parton) structure of nucleons and nu¬ 
clei. Structure functions extracted from inclusive deep- 
inelastic scattering (DIS) experiments display the central 
features of quantum chromodynamics (QCD) — asymp¬ 
totic freedom at short distances (via structure function 
scaling and its violation) and confinement at large dis¬ 
tance scales (via the momentum dependence of parton 
distributions). 

Since the late 1960s, DIS experiments have yielded an 
impressive data set that maps nucleon structure functions 
over several orders of magnitude in the Bjorken scaling 
variable, x, and the squared four-momentum transfer, 
Q 2 . These data, supplemented by cross sections from 
hadronic collisions and other high-energy processes, have 
enabled a detailed picture of the parton distribution func¬ 
tions (PDFs) of the nucleon to emerge through global 
QCD analyses (see Refs. mu and references therein). 

At lower energies, where nonperturbative quark-gluon 
interactions are important and the inclusive lepton- 
nucleon cross section is dominated by nucleon resonances, 
the structure functions reveal another intriguing feature 
of QCD, namely, quark-hadron duality. Here, the low 
energy cross section, when averaged over appropriate en¬ 
ergy intervals, is found to resemble the high energy result, 
whose Q 2 dependence is described by perturbative QCD. 
In this context, quark-hadron duality provides a unique 
perspective on the relationship between confinement and 
asymptotic freedom, and establishes a critical link be¬ 
tween the perturbative and nonperturbative regimes of 
QCD. 

In the framework of QCD, quark-hadron duality can 
be formally interpreted in terms of structure function mo¬ 


ments [3j. From the operator product expansion (OPE), 
the moments can be expressed as a series in 1 /Q 2 , with 
coefficients given by matrix elements of local quark-gluon 
operators of a given twist [3J. The leading (twist 2) term 
corresponds to scattering from a single parton, while 
higher twist terms correspond to multi-quark and quark- 
gluon interactions. At low Q 2 the resonance region makes 
a significant contribution to the structure function mo¬ 
ments, so that here one might expect a strong Q 2 depen¬ 
dence of the moments arising from the higher twist terms 
in the OPE. In practice, however, the similarity of the 
structure function moments at low Q 2 and the moments 
extracted from high energy cross sections suggests the 
dominance of the leading twist contribution. The com¬ 
bined higher twist, multi-parton contributions appear to 
play a relatively minor role down to scales of the order 
Q 2 ~ 1 GeV 2 . 

This nontrivial relationship between the low-energy 
cross section and its deep-inelastic counterpart was first 
observed by Bloom and Gilman in the early DIS 

measurements that were instrumental in establishing 
structure function scaling. More recently, the availability 
of extensive, precise structure function data from Jeffer¬ 
son Lab and elsewhere, over a wide range of kinemat¬ 
ics, has opened up the possibility for in-depth studies of 
quark-hadron duality. Duality has now been observed 
in the proton F 2 and Fj, structure functions [M2], the 
F 2 structure function of nuclei [13], the spin-dependent 
gi structure functions of the proton and 3 He P3HI5], the 
individual helicity-1/2 and 3/2 virtual photoproduction 
cross sections for the proton mi and in parity-violating 
electron-deuteron scattering [151 . 

To establish the dynamical origin of quark-hadron du¬ 
ality in the nucleon requires one to also study the low-Q 2 
structure of the neutron. Models based on four-quark 
higher twist contributions to DIS suggest that duality 
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in the proton could arise from accidental cancellations 
between quark charges, which would not occur for the 
neutron nsi- Unfortunately, the absence of high-density 
free neutron targets means that essentially all informa¬ 
tion on the structure functions of the neutron has had 
to be derived from measurements on deuterium. Typi¬ 
cally, the deuterium data are corrected for Fermi smear¬ 
ing and other nuclear effects [2fiH2()| . which introduces 
an element of model dependence into the extraction pro¬ 
cedure. This is particularly problematic in the nucleon 
resonance region, where Fermi motion effects leads to 
significant smearing of the resonant structures. The ex¬ 
istence of duality in the neutron F 2 structure function 
was suggested recently [20j in an analysis which used an 
iterative deconvolution method m to extract neutron 
resonance spectra from inclusive proton and deuteron F 2 
data [9]. A model independent confirmation of duality in 
the neutron, however, was to date not possible. 

Recently, a new experimental technique, based on 
spectator nucleon tagging [25], has been used to extract 
the free neutron F 2 structure function [29]. By detecting 
low-momentum protons at backward angles in electron 
deuteron scattering, the BONuS experiment at Jeffer¬ 
son Lab measured F 2 for the neutron in both the reso¬ 
nance and DIS regions, with minimal uncertainty from 
nuclear smearing and rescattering corrections [30]. In 
the present work, we use the BONuS data to quantita¬ 
tively measure for the first time the degree to which du¬ 
ality holds for the F 2 structure function of the free neu¬ 
tron. Because the results reported here use data from 
an experimentally-isolated neutron target, one expects 
significantly reduced systematic uncertainties compared 
with those in the model-dependent analysis of inclusive 
deuterium data [20] . 

For the theoretical analysis of duality we use the 
method of truncated structure function moments [H], 
which were applied to the resonance region F 2 proton 
data by Psaker et al. [35]. Here, the n-th truncated 
moment of the F 2 structure function is defined as 


MN(x m i n , x max , Q ) — 


“ J-min 


dxx"- 2 F 2 {x,Q 2 ), (1) 


where the integration over x is restricted to an interval 
between ir m i n and x max . This method avoids extrapola¬ 
tion of the integrand into poorly mapped kinematic re¬ 
gions, and is particularly suited for the study of duality 
where an x interval can be defined by a resonance width 
around an invariant mass W 2 = M 2 +Q 2 (l — x)/x, where 
M is the nucleon mass. As the position of the resonance 
peak varies with x for different Q 2 values, the values for 
;r m i n and x max evolve to the appropriate invariant mass 
squared region. For the BONuS data, we consider four 
ranges in W 2 , corresponding to the three prominent res¬ 


onance regions as well as the combined resonance region, 

1.3 < W 2 <1.9 GeV 2 [1st (or A) region], 

1.9 < W 2 < 2.5 GeV 2 [2nd region], 

2.5 < W 2 < 3.1 GeV 2 [3rd region], (2) 

1.3 < W 2 < 4.0 GeV 2 [total resonance]. 

After reviewing the BONuS experiment in Sec. [TT] the 
results for several low truncated moments (correspond¬ 
ing to N = 2, 4 and 6) of the neutron F 2 structure 
function are presented in Sec. m The implications of 
the new data for local quark-hadron duality and its vi¬ 
olation are discussed by comparing with recent global 
PDF parametrizations and previous model-dependent 


data analyses (Sec. Ill AI. The isospin dependence of lo¬ 
cal duality is studied by comparing the neutron moments 
with corresponding moments of the proton F 2 structure 
function (Sec. Ill B). Finally, we summarize our results 
and conclusions in Sec. hyi 


II. THE BONUS EXPERIMENT 

The results reported here rely on a novel experimen¬ 
tal technique aimed at eliminating or substantially re¬ 
ducing the theoretical uncertainties involved in extract¬ 
ing neutron data from nuclear targets. The BONuS 
(Barely Off-shell Nucleon Structure) experiment at Jef¬ 
ferson Lab [281430] used a Radial Time Projection Cham¬ 
ber (RTPC) to detect low momentum spectator protons 
produced in electron-deuterium scattering in conjunction 
with electrons detected with CLAS [SO] in Hall B. The 
tagging technique essentially eliminates effects of Fermi 
smearing 133, while the restriction to backward low- 
momentum spectator protons minimizes final state in¬ 
teractions [381440] and ensures that the neutron is nearly 
on-shell [29j [41] . 

The BONuS experiment ran in 2005 and acquired 
electron-deuteron scattering data at two electron beam 
energies, E = 4.223 and 5.262 GeV. The RTPC consisted 
of three layers of gas electron multipliers surrounding a 
thin, pressurized gas deuterium target, and was able to 
detect protons with momenta as low as 70 MeV. The 
experiment and data analysis are described in detail in 
Ref. [30]. Ratios of neutron to proton F 2 structure func¬ 
tions and the absolute neutron F 2 structure function were 
extracted over a wide kinematic range and for spectator 
proton momenta between 70 and 100 MeV. The total 
systematic uncertainty in the neutron structure function 
extracted was 8.7% [30], with an overall 10% scale uncer¬ 
tainty. 

The kinematic coverage, shown in Fig. |T] (with the 
4.223 and 5.262 GeV data combined), extends from the 
threshold to the deep-inelastic region. The curves in 
Fig. 0 represent the fixed-IV 2 thresholds for the four 
mass regions considered. Typical neutron F 2 spectra 
are shown in Fig. 1 for Q 2 = 1.2 and 2.4 GeV 2 , with the 
data restricted to spectator proton angles greater than 
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FIG. 1: Kinematic coverage of the BONuS data. The solid 
lines denote the fixed-IU 2 thresholds for the four final state 
mass regions in Eq. ([ 2 ]), from W 2 = 1.3 to 4.0 GeV 2 . 



FIG. 3: Second (N = 2) neutron truncated moments 
versus Q 2 for the four resonance regions in Eq. (§ [la¬ 
beled as “first”, “second”, “third” and “total”]. The mo¬ 
ments obtained from the BONuS data (filled circles) are com¬ 
pared with moments computed from the ABKM global PDF 
parametrization [42] including target mass and higher twist 
corrections (shaded rectangles). 


III. TRUNCATED MOMENTS AND LOCAL 
QUARK-HADRON DUALITY 

Because the kinematic variables Q 2 , x and W 2 are not 
independent, a range in W 2 at fixed Q 2 implies a corre¬ 
sponding range in x. This makes possible a straightfor¬ 
ward integration of the experimental F,™ structure func¬ 
tion data to obtain the truncated moments M n in Eq. ([l]) . 
To minimize the model dependence, we evaluate the in¬ 
tegrals based solely on the experimentally measured data 
points, without using any interpolating or extrapolating 
function. 


FIG. 2: Representative neutron F 2 "' structure function spec¬ 
tra from the BONuS experiment [30] at Q 2 = 1.2 GeV 2 
(top panel) and Q 2 = 2.4 GeV 2 (bottom panel). The open 
(filled) circles represent data for a beam energy of E = 4.223 
(5.262) GeV. The solid curve is computed from the ABKM 
global PDF parametrization [42] including higher twist effects 
and target mass corrections. 


100° relative to the momentum transfer, and proton mo¬ 
menta between 70 and 100 MeV. The BONuS results are 
compared with the ABKM global fit [12] to deep-inelastic 
and other high-energy scattering data, with the inclusion 
of target mass corrections and higher twist effects. The 
qualitative agreement between the parametrization and 
data suggests evidence for quark-hadron duality, which 
we explore more quantitatively in the following sections. 


A. Truncated neutron moments 

The second (N = 2) truncated neutron moments, Mg, 
obtained from the BONuS data are shown in Fig. [3] as 
a function of Q 2 for the four W 2 intervals defined in 
Eq. ([ 2 ]). The numerical values for the moments are also 
listed in Table [I] The quoted errors take into account 
the experimental statistical and systematic uncertainties 
added in quadrature, but do not the include the 10% 
scale uncertainty. 

The experimental moments are compared with the 
moments calculated from the ABKM global PDF 
parametrization [42] , including finite-Q 2 corrections from 
the target mass and an x-dependent parameterization of 
the leading ((D(l/Q 2 )) higher twist effects. The latter are 
needed in order to obtain a more quantitative descrip¬ 
tion of duality in the low-Q 2 region, to which the struc¬ 
ture functions from the global fits (which are primarily 
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Ah [xio~ 3 ] 


Q 2 [GeV"] 

1st 

2nd 

3rd 

total 

1.0 

31.5 ± 1.1 

16.4 ± 0.4 

14.1 

± 0.3 

76.7 

± 

1.2 

1.2 

23.5 ± 0.5 

15.3 ± 0.3 

13.0 

± 0.3 

67.4 

± 

0.6 

1.4 

17.7 ± 0.4 

13.5 ± 0.2 

12.1 

± 0.3 

57.7 

± 

0.5 

1.7 

12.3 ± 0.3 

10.7 ± 0.2 

10.5 

± 0.2 

46.7 

± 

0.5 

2.0 

8.4 ± 0.2 

8.5 ± 0.2 

9.0 

± 0.2 

38.4 

± 

0.4 

2.4 

5.8 ± 0.2 

6.1 ± 0.1 

7.1 

± 0.3 

29.5 

± 

0.4 

2.9 

3.4 ± 0.1 

4.5 ± 0.1 

5.3 

± 0.3 

21.5 

± 

0.4 

3.4 

2.1 ± 0.1 

3.1 ± 0.1 

4.0 

± 0.2 

15.8 

± 

0.3 

4.1 

1.3 ± 0.1 

1.7 ± 0.1 

2.8 

± 0.1 





TABLE I: Second (N = 2) truncated moments (in units of 
10~ 3 ) of the neutron 7*2 structure function from the BONuS 
data for the W 2 regions in Eq. |2|. The errors are a quadra¬ 
ture sum of statistical and systematic uncertainties, but do 
not include the overall 10% normalization uncertainty. 



constrained by high energy data) are extrapolated. The 
comparison shows generally very good agreement in the 
second and third resonance regions, and in the total in¬ 
tegrated W 2 interval, while the ABKM results underes¬ 
timate the data somewhat in the A resonance region. 

The corresponding higher order truncated moments 
(for N = 4 and N = 6) are listed in Tables [II] and III re¬ 
spectively. Comparison with the ABKM fit(not shown) 
reveals a similar pattern as for the N = 2 moments, al¬ 
though the deviation in the lowest-VF interval is more 
pronounced, especially at low Q 2 , because of the greater 
weighting given to the high-a; region in the higher mo¬ 
ments. 


Ah [x 10~ :i ] 


Q 2 [GeV"] 

1st 

2nd 

3rd 

total 

1.0 

11.58 

it 

0.43 

3.09 

± 

0.08 

1.69 

± 

0.04 

17.49 

± 

0.44 

1.2 

9.80 

± 

0.21 

3.51 

± 

0.06 

1.95 

± 

0.04 

16.78 

± 

0.22 

1.4 

8.11 

± 

0.17 

3.60 

± 

0.06 

2.17 

± 

0.04 

15.61 

± 

0.19 

1.7 

6.27 

± 

0.14 

3.40 

± 

0.06 

2.33 

± 

0.05 

14.01 

± 

0.17 

2.0 

4.67 

± 

0.14 

3.08 

± 

0.06 

2.36 

± 

0.06 

12.45 

± 

0.17 

2.4 

3.48 

± 

0.11 

2.54 

± 

0.06 

2.20 

± 

0.08 

10.59 

± 

0.15 

2.9 

2.22 

± 

0.10 

2.11 

± 

0.07 

1.93 

± 

0.09 

8.52 

± 

0.16 

3.4 

1.44 

± 

0.09 

1.58 

± 

0.07 

1.64 

± 

0.08 

6.72 

± 

0.15 

4.1 

0.95 

± 

0.08 

0.98 

± 

0.07 

1.29 

± 

0.06 





TABLE II: As in Table [I] but for the N = 4 moment. 


Note that while early phenomenological analyses of 


A/s [x!0~ 3 ] 


Q 2 [GeV"] 

1st 

2nd 

3rd 

total 

1.0 

4.39 

± 

0.18 

0.60 

± 

0.01 

0.20 

± 

0.01 

5.28 

± 

0.18 

1.2 

4.19 

± 

0.10 

0.82 

± 

0.01 

0.30 

± 

0.01 

5.45 

± 

0.10 

1.4 

3.79 

± 

0.09 

0.98 

± 

0.02 

0.39 

± 

0.01 

5.38 

± 

0.09 

1.7 

3.24 

± 

0.08 

1.09 

± 

0.02 

0.52 

± 

0.01 

5.17 

± 

0.08 

2.0 

2.62 

± 

0.08 

1.13 

± 

0.02 

0.62 

± 

0.02 

4.82 

± 

0.09 

2.4 

2.12 

± 

0.07 

1.06 

± 

0.02 

0.68 

± 

0.02 

4.41 

± 

0.08 

2.9 

1.45 

± 

0.07 

1.00 

± 

0.03 

0.71 

± 

0.03 

3.77 

± 

0.08 

3.4 

0.99 

± 

0.07 

0.82 

± 

0.04 

0.67 

± 

0.03 

3.14 

± 

0.09 

4.1 

0.68 

± 

0.06 

0.56 

± 

0.04 

0.60 

± 

0.03 


— 



TABLE III: As in Table [I] but for the N = 6 moment. 


FIG. 4: Ratios of truncated moments of the neutron F 2 struc¬ 
ture function from the BONuS data to those computed from 
the ABKM global PDF parametrization [42] including finite- 
Q 2 effects (filled circles) as a function of Q 2 for the four W 2 
intervals in Eq. §. The empirical moments are compared 
with the results of the model-dependent analysis of inclusive 
DIS data (201 (open circles), and with ratios computed from 
the CJ12 distributions 33], with leading twist only (dotted 
lines) and including hnite-Q 2 effects (solid lines). All ratios 
are taken relative to the ABKM moments. 


quark-hadron duality typically compared resonance re¬ 
gion data at low Q 2 with scaling functions extrapolated 
from fits to high-IK cross sections 0 0 , more recent 
quantitative analyses 0HD] have emphasized the need to 
take into account the Q 2 dependence in the high-kF data, 
including both leading and higher twist contributions. 
This is especially important in the high-x region, where 
the separation between the leading and higher twists is 
more model dependent due to the absence of high-Q 2 
measurements, and comparison of resonance region data 
with the total extrapolated structure functions reveals an 
enhanced persistence of duality down to lower values of 

Q 2 . 

To study local quark-hadron duality in detail, we form 
ratios of the truncated moments of obtained from the 
BONuS data to the moments computed from the ABKM 
reference structure function [42], over the same range of 
x. The ratios for the M% moments are shown in Fig. [4] 
as a function of Q 2 for the four invariant mass regions 
in Eq. ([ 2 ]). The ratios for the second, third and total 
resonance regions are close to unity, to within ~ 10% 
over nearly the entire range of Q 2 = 1 — 4 GeV 2 , and ex¬ 
hibit weak scale dependence. This points to a dramatic 
confirmation of the validity of local duality for the neu¬ 
tron in these regions. In the first resonance region, the 
A resonance is ~ 20% — 30% larger than the PDF-based 
fit, but still displays a similar Q 2 behavior. This could 
be interpreted as either a stronger violation of local du¬ 
ality in the A region, which may be expected at lower 
W, or possibly underestimated strength of the ABKM 
parametrization in the large-x regime, to which this W 
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region corresponds. While this is difficult to disentan¬ 
gle experimentally, duality is expected to hold to better 
accuracy for integrals obtained over regions containing 
multiple final states. 

The confirmation of the approximate validity of dual¬ 
ity in F£ from the BONuS data disfavors the suggestion 
m that duality occurs in the proton because of acciden¬ 
tal cancellations of quark charges associated with higher 
twist, four-quark operators, and disagrees with the pre¬ 
diction that duality should therefore not be seen in the 
neutron. This conclusion was also reached in the model- 
dependent analysis by Malace et al. [20], who studied 
duality in the neutron by extracting the F£ structure 
function from inclusive DIS data using phenomenologi¬ 
cal deuteron wave functions and an iterative deconvolu¬ 
tion procedure EZj. Overall, the BONuS data are in good 
agreement with the earlier results [201 , within the experi¬ 
mental uncertainties, although they appear to lie system¬ 
atically higher in the A region. This may be associated 
with the nuclear corrections in the deuteron, which are 
subject to greater uncertainties at the largest x (small¬ 
est W) values, or a systematic bias of the subtraction 
method in relation to the various theoretical assumptions 
and models [21] , 

The relevance of large-a; uncertainties and Hnite-Q 2 
corrections in global PDF fits is illustrated in Fig. [4j 
where the experimental and computed ABKM moments 
are also compared with the moments calculated from the 
CTEQ-Jefferson Lab (CJ) global PDF parametrization 
[33] with and without higher twist corrections. While the 
ratio of the ABKM and C J moments is close to unity over 
the entire range of Q 2 considered when finite-Q 2 effects 
are included, the deviation from unity of the ratio com¬ 
puted from only the leading twist components of the CJ 
fit can be up to 30% — 40% for the integrated resonance 
region, and up to twice as much for the A region. This 
suggests an important role played by the finite-Q 2 cor¬ 
rections to the scaling functions in effecting the cancella¬ 
tions between the individual resonance regions necessary 
for the realization of quark-hadron duality pHT - UH] . 

However, even incorporating finite-Q 2 corrections, 
global PDF fits can differ significantly in the large-x (low- 
W) regime. Because of the Q 2 and W 2 cuts applied to 
the global data set, PDFs in the largest-a; regions relevant 
for this analysis are essentially unconstrained, resulting 
in large uncertainties in the extrapolated x —>■ 1 behavior 

@3- 


B. Isospin dependence 

The stronger violation of local duality in the A re¬ 
gion is also evident in the ratio of neutron to proton 
truncated moments, shown in Fig. [5] compared with the 
reference ABKM parametrization [42] that was used to 
compute both the proton and neutron moments. To ob¬ 
tain the empirical proton truncated moments in the res¬ 
onance region, the Christy-Bosted global fit [IS] was uti- 



FIG. 5: Ratios M£/M% of neutron to proton truncated mo¬ 
ments of the F 2 structure function versus Q 2 , for the four W 
regions in Eq. <§■ The BONuS results (filled circles) are com¬ 
pared with the moments computed from the ABKM global 
PDF parametrization including target mass and higher twist 
(solid lines) corrections. In both cases the proton moments 
are evaluated from the same ABKM fit m- 

lized. (Duality in the proton structure function moments 
themselves was studied in detail in previous analyses [5] , 
and generally confirmed at the 10% — 15% level for the 
N = 2 moment when integrated over the entire resonance 
region.) 

The significant duality violation in the neutron/proton 
ratio observed in the A region can be understood from 
the isovector nature of the A-isobar and the relatively 
small nonresonant background on which it sits. In the 
limit of exact isospin symmetry, the transitions from a 
ground state nucleon to an isospin-3/2 resonance would 
be identical for protons and neutrons. Nonresonant back¬ 
ground and isospin symmetry breaking contributions give 
rise to differences between proton and neutron moments, 
but these are typically very small in the A region. In 
contrast, the proton and neutron deep-inelastic struc¬ 
ture functions (either leading twist only or with higher 
twist corrections) in the A region are expected to be 
quite different, since at large x the neutron structure 
function is strongly suppressed relative to the proton, 
F -2 <C F -2 [49j 50]. The fact that the experimental 
M 2 /ratio in the high-x region lies somewhat higher 
than the DIS parametrization (even more pronounced 
than in Ref. m) is therefore consistent with these ex¬ 
pectations. 

A similar comparison of the neutron to proton mo¬ 
ments in the second and third resonance regions in 
Fig. [5] shows significantly better agreement with the DIS 
parametrization. Based on simple quark models and as¬ 
suming magnetic coupling dominance, one would expect 
the resonance contribution to the neutron moments to 
underestimate the DIS moment in the second resonance 
region. This is due to the small couplings to octet spin- 
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1/2 states. In contrast, according to Refs. [1SJ HU [FT] the 
proton moments would overestimate the DIS results in 
the second and third W intervals in Eq. §■ While there 
was some evidence for such a pattern from the earlier, 
model-dependent analysis of inclusive data m, there is 
no indication from the BONuS results of a suppression 
in the second resonance region. The slightly larger over¬ 
all magnitude of the neutron moments compared with 
Ref. m brings the present results into excellent agree¬ 
ment with the DIS moments in the second region, with a 
small enhancement in the third region. The correspond¬ 
ing enhancement of the proton data in the third reso¬ 
nance region relative to the ABKM fit El then results 
in essentially no deviation of the neutron to proton ratio 
here, as illustrated in Fig. [5] Finally, for the total inte¬ 
grated region between threshold and W = 2 GeV, the 
empirical M%/M% ratio is slightly above the DIS result 
mostly because of the large enhancement of the data in 
the A region. 

IV. CONCLUSION 

In this work we have investigated local quark hadron 
duality in the neutron structure function based on data 
from the BONuS experiment at Jefferson Lab [29] 30] . 
which used a novel experimental technique to create an 
effective neutron target by tagging low momentum spec¬ 
tator protons in electron-deuterium scattering. The spec¬ 
tator tagging technique provides smaller systematic un¬ 
certainties compared with the traditional method of sub¬ 
tracting smeared hydrogen data and from inclusive deu¬ 
terium structure functions, using model assumptions for 
the nuclear corrections. 

We have evaluated the N = 2, 4 and 6 truncated mo¬ 
ments of the neutron F£ structure function for the three 
standard nucleon resonance regions and the total inte¬ 
grated resonance region up to W = 2 GeV, over the 
range Q 1 2 3 4 = 1.0 to 4.1 GeV 2 . Comparison of the exper¬ 
imental moments with moments computed from global 
parametrizations of PDFs fitted to deep-inelastic and 
other high energy scattering data, as well as with the 
corresponding truncated moments for the proton, reveals 
a dramatic confirmation of local duality for the neu¬ 
tron in the second, third and total resonance regions to 
better than 10% for the lowest moment. The stronger 
(~ 20% — 30%) violation of duality in the A region is con¬ 
sistent with the expectations based on isospin symmetry 


for the isovector transition amplitudes and the behavior 
of the F% /F% ratio at large x [43] [50]. 

The confirmation of local duality in the neutron disfa¬ 
vors the model m in which duality in the proton arises 
through accidental cancellations of quark charges asso¬ 
ciated with higher twist, four-quark operators, which 
would predict strong duality violations in the neutron. 
Rather, it suggests a dynamical origin of duality in which 
cancellations among nucleon resonances produce a higher 
degree of duality over the entire resonance region, with 
stronger violations locally |T51 [151 [5T| . On the other 
hand, detailed comparisons between the empirical trun¬ 
cated moments and DIS parametrizations in the individ¬ 
ual resonance regions suggest a pattern of duality viola¬ 
tion that is more involved than that predicted by simple 
spin-flavor symmetric quark models with magnetic cou¬ 
pling dominance. 

Our results also confirm and refine the findings of ear¬ 
lier model-dependent studies [20] of duality in the neu¬ 
tron in which the neutron structure was extracted from 
inclusive proton and deuteron data assuming a model 
for the nuclear corrections and an iterative deconvolu¬ 
tion procedure m- In particular, the BONuS moments 
are found to lie slightly higher than the earlier results, 
especially in the A region, but with a similar Q 2 depen¬ 
dence. 

In the future, the spectator tagging technique will be 
used at Jefferson Lab with an 11 GeV electron beam to 
extend the kinematical coverage of K/ measurements to 
higher values of x and Q 2 [52;. As well as providing 
more stringent constraints on the leading twist PDFs in 
the limit x — » 1, the new data will allow more definitive 
tests of local quark-hadron duality for the neutron over 
a greater range of Q 2 . 
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